The Harris-Luck criterion for random lattices 
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The Harris-Luck criterion judges tlie relevance of (potentially) spatially correlated, quenched 
disorder induced by, e.g., random bonds, randomly diluted sites or a quasi-periodicity of the lattice, 
for altering the critical behavior of a coupled matter system. We investigate the applicability of 
this type of criterion to the case of spin variables coupled to random lattices. Their aptitude to 
alter critical behavior depends on the degree of spatial correlations present, which is quantified 
by a wandering exponent. We consider the cases of Poissonian random graphs resulting from the 
Voronoi'-Delaunay construction and of planar, "fat" 4>'^ Feynman diagrams and precisely determine 
their wandering exponents. The resulting predictions are compared to various exact and numerical 
results for the Potts model coupled to these quenched ensembles of random graphs. 

PACS numbers: 75.10.Hk, 75.40.Mg, 75.50.Lk 



I. INTRODUCTION 

The concept of quenched disorder coupling to the local 
energy density has evolved to be the main framework for 
the modeling of the types of randomness found in real 
physical systems^'^. In the presence of frustration this 
includes the vast field of spin glasses, which in the past 
decades has attracted an enormous amount of analytical 
and numerical research, see, e.g., Refs. 2-5. Here, we 
will be concerned with the simpler case of models with 
purely ferromagnetic couplings. The first investigations 
of this problem have considered defects distributed in the 
system completely at random^'^, realized in lattice mod- 
els, e.g., as random variation of the coupling strengths or 
random deletion of bonds or lattice sites*"^°. 

The relevance of this random-bond or dilution type of 
disorder for the universal behavior of spin systems has 
been the subject of much research*'^^"^^. For the case 
of models undergoing a continuous phase transition on 
regular lattices, Harris* argued that for models with a 
specific-heat exponent a < the fluctuations in the local 
transition temperature induced by the disorder degrees 
of freedom are not strong enough to alter universal fea- 
tures of the model such as the critical exponents. Albeit 
not originally claimed by Harris*, for the converse case of 
a positive specific-heat exponent a a significant change 
of the system's behavior was expected. The precise ef- 
fect of such a relevant perturbation has been the subject 
of some debate^'*. While it was originally believed that 
the transition temperature fluctuations might smoothen 
out the phase transition so far as to completely destroy 
it, it was later on realized that this, in fact, does not 
happen and, instead, the system experiences a cross-over 
from the pure fixed point to a new, disorder fixed point, 
resulting in a new set of critical exponents and further 
universal properties such as amplitude ratios^'^^'^'''^*"^°. 
This scenario has been especially thoroughly analyzed for 
the case of the q ~ 2,3,4 Potts models in two dimensions, 
where results from perturbative mcthods"'^'^"^^ agree well 
with the outcome of numerical simulation studies^°'^^"^'^. 
From this general observation of a smoothening effect of 



disorder on phase transitions, one might expect that for 
systems exhibiting a first-order transition in the regu- 
lar case, disorder of the random-bond type might soften 
the transition to a continuous one^^. For the case of 
two dimensions it could be rigorously established that 
even an infinitesimally small amount of disorder suffices 
to indeed induce this behavior^"*^^^; in three dimensions, 
sufficiently strong disorder coupling to the local energy 
density is numerically found to soften first-order phase 
transitions to second-order ones, see, e.g., Refs. 28-30. 

Obviously, for many physical systems the assumption 
of an uncorrelated, isotropic distribution of defects is not 
an adequate description. Instead, due to various rea- 
sons the distribution of defects is spatially correlated. 
This effect occurs isotropically due to long-range inter- 
actions between the non-magnetic ions, or in the form 
of line or higher-dimensional defects^ In these cases, 
the reasoning leading to Harris' relevance criterion is no 
longer directly applicable, but can be generalized accord- 
ingly. For the case of algebraically decaying correlations 
one finds a relevance threshold depending on the dimen- 
sion of the defects as well as the power of the decay of 
the correlations^^''^^. For a different model, not covered 
by the random-bond paradigm, namely the co-ordination 
number non-periodicity found in quasi-crystals and other 
aperiodic structures, Luck^^ formulated a relevance cri- 
terion that includes the situations discussed as special 
cases. There, the "break-even point" for the relevance of 
randomness in terms of the specific-heat exponent a is 
shifted from its uncorrelated value ac = to somewhere 
in the region — oo < etc < 1, depending on the strength 
of spatial correlations of the disorder degrees of freedom 
measured by a geometrical fluctuation or wandering expo- 
nent. A structure conceptually related to these aperiodic 
models is given by different varieties of random graphs 
such as Poissonian Voronoi-Delaunay triangulations'^^''^^ 
or the planar, combinatorial triangulations encountered 
in the dynamical triangulations approach to quantum 
gravity resp. the dual planar and orientable (f>^ Feynman 
diagrams'^*. Although spin models on quenched ensem- 
bles of these graphs have been considered in numerical 



FIG. 1: Patch of a Poissonian Voronoi-Delaunay triangulation 
of the plane. The spherical, shaded dots denote the randomly 
distributed generators. The solid lines enclosing these gener- 
ators define convex polygonal areas whose points are closer to 
the origin than to any other generator. They form the network 
of Voronoi cells, which results in a three-valent graph whose 
vertices are depicted by the dark boxes. Its geometrical dual, 
consisting of the generators and the connecting dashed lines, 
is the Delaunay triangulation. 

simulation studies'^^"^^, no connection has been made as 
yet with the predictions of Luck^^ for general systems 
with connectivity disorder. In particular, the geomet- 
rical fluctuation exponents appearing in this relevance 
criterion have not been determined for the case of these 
random graphs. 

The rest of this paper is organized as follows. In Sec- 
tion II we introduce the ensembles of graphs to be consid- 
ered, elaborate on their generation in a computer exper- 
iment and review some of their known properties. Sec- 
tion III is devoted to a presentation of a version of the 
Harris-Luck relevance criterion suitable to be applied to 
the case of random graphs considered here and a discus- 
sion of its connection with previous results. In Section 
IV we introduce different methods to precisely determine 
the wandering exponent for both graph types and present 
the results obtained. Finally, Section V contains our con- 
clusions. 



II. RANDOM GRAPHS 

In the following we present construction techniques 
and properties of two different kinds of two-dimensional 
topological graphs whose randomness is solely encoded 
in the degree distribution of their vertices, resulting in 
a topological connectivity disorder. In contrast to the 
generic random graphs discussed in the context of scale- 
free networks, small- world models, etc.'^'^'^^, the ensem- 
bles of random structures considered here are not fully 
determined by their degree or co-ordination number dis- 
tributions, but additionally exhibit a well-defined topol- 
ogy and long-range correlations of their disorder degrees 



FIG. 2: (a) The link-flip move on two adjacent triangles of a 
dynamical triangulation (dashed lines). The solid lines denote 
the corresponding dual, three-valent cj)^ graphs, (b) Example 
of a spherical dynamical triangulation embedded into three- 
dimensional space. Note that the requirement of equilater- 
ality of the triangles is only approximately fulfilled for the 
embedding shown, since there is no generic exact embedding 
algorithm for performing it. 



of freedom. 



A. Poissonian Voronoi'-Delaunay Triangulations 

Irregular cell structures or froths appear in a large 
variety of natural systems, such as foams or biological 
tissues'^^. The common inverse problem of constructing 
such a connectivity or cell structure from a given irregular 
arrangement of vertices (so-called generators) is solved by 
the Voronoi-Delaunay construction'^^. In two dimensions, 
the Voronoi cell of a given generator is a convex polygon 
around it, enclosing the part of its neighborhood which 
is closer to it than to any other generator, cf. Fig. 1. 
This is in complete analogy with the notion of a Wigner- 
Seitz elementary cell in crystallography. This construc- 
tion results in the three-valent Voronoi graph and the 
dual Delaunay triangulation. If the generators are lo- 
cated completely at random, as is the case in the exam- 
ple of Fig. 1 , the resulting graph is known as Poissonian 
Voronoi-Delaunay triangulation"^^ . 

These lattices are randomly disordered in several re- 
spects: edge lengths, cell volumes etc. vary, as well as 
the co-ordination numbers q of the vertices of the De- 
launay triangulation. To facilitate comparison with the 
second type of random graphs to be introduced below, we 
here restrict ourselves to the latter aspect of variable co- 
ordination numbers, i.e., we do not take any effects from 
variable lengths or areas into account. Thus, the distri- 
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bution P{q) of co-ordination numbers is the only random 
variable involved. Additionally, to eliminate surface ef- 
fects, the generators are randomly placed on the surface 
of a sphere instead of a patch of the plane. Thus, the 
resulting random graphs are triangulations of spherical 
topology. 

From the Euler relations, the average co-ordination 
number is a topological invariant for a fixed number of 
triangles in two dimensions, given by'^* 

^"-^E'^'^e^, (1) 

i 

for any closed, spherical triangulation, where N denotes 
the number of triangles. In the limit of infinite triangu- 
lations, N — > cxD, one thus obviously has (q) = (qi) — 6. 
The second moment of q is not exactly known, but is 
numerically found to be"^^'^^ 

/i2 = {qf) - « 1.781, (2) 

as ^ CXD. It turns out that the random variables qi are 
not independently distributed, but are reflecting a spatial 
correlation of the disorder degrees of freedom in addition 
to the trivial correlation induced by the constraint (1). 
The form of these correlations for nearest-neighbor ver- 
tices is commonly described by the Aboav-Weaire law"^^, 
which states that the total expected number of edges of 
the neighbors of a g-sided cell, qm{q), should vary lin- 
early with q, 

qm{q) ^ {6 - a)q + b, (3) 

where a and b are some parameters. In turns out, how- 
ever, that Eq. (3), albeit being a good effective descrip- 
tion for a large variety of cell systems including the 
case of Poissonian Voronoi-Delaunay triangulations, is 
merely a leading-order and not an exact property of these 
systems^''. 

B. Dynamical Triangulations 

An ensemble of planar triangulations with proper- 
ties very different from those of the Poissonian Voronoi- 
Delaunay lattices is given by the so-called dynamical tri- 
angulations model, which has been used as a construc- 
tive approach to discrete Euclidean quantum gravity in 
two dimensions'^. This ensemble of combinatorial tri- 
angulations is defined as that of all gluings of equilat- 
eral triangles to closed surfaces of a given topology (such 
as, e.g., that of a sphere), where all gluings are counted 
with equal probability. A dynamical way of generating 
such graphs is the repeated application of so-called link- 
flip moves which re- wire a given network [see Fig. 2(a)]. 
It is known that this procedure converges to a limit- 
ing distribution with specific properties described below. 
Since this flip dynamics introduces large temporal cor- 
relations between subsequent random triangulations, we 




FIG. 3: Comparison of the degree distributions P{q) of Pois- 
sonian Delaunay triangulations and dynamical triangulations 
in the limit of an infinite number of triangles. The results are 
taken from Refs. 49,50. 

do not actually use it to generate the graph instances 
needed below, but instead revert to a recursive inser- 
tion technique known to yield the same graph distribu- 
tion, but independent graph realizations in each step^^. 
In this construction, no edge length differences are in- 
volved and, again, the randomness is solely encoded in 
the degree distribution P{q) of the vertices. Although 
these objects are entirely defined in terms of their intrin- 
sic connectivity properties, an embedding of an exam- 
ple dynamical triangulation into three-dimensional space 
is shown for illustration purposes in Fig. 2(b). Techni- 
cally, the geometrical duals of these triangulations are 
given by the ensemble of planar, "fat" (i.e., orientable) 
(j)'^ Feynman diagrams without tadpoles and self-energy 
insertions. The statistics of these objects can be explic- 
itly treated by means of matrix models^^'^'^, leading to 
exact solutions for the co-ordination number distribution 
P{q) and many other properties of the modeP^'*^". Com- 
pared to the Poissonian Voronoi-Delaunay model, fluctu- 
ations are much more pronounced for these graphs, and it 
can be shown that the variance of co-ordination numbers 
approaches'"'^"* 

/i2 = 21/2 (4) 

'AS N oo, whereas Eq. (1) still holds. Also, the Aboav- 
Weaire law (3) correctly describes the leading-order be- 
havior of the nearest-neighbor correlations^'*. 

Considering more geometric properties of the graphs, 
the differences between both ensembles become even 
more pronounced. Especially, the quantum gravity 
graphs of the dynamical triangulations model can be 
shown to be highly fractal, being self-similarly composed 
of "baby universes" branching off from the main sur- 
face, i.e., macroscopical subgraphs attached to the main 
body by only a few links'' . This fractal structure leads 
to an exceptionally large internal Hausdorff dimension 
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of dh = 4: as compared to the topological dimension of 
two^^'''''', whereas for the Poissonian Voronoi-Delaunay 
triangulations the Hausdorff dimension remains at the 
trivial value of the topological dimension. The differ- 
ences between the two graph ensembles are also very 
strikingly seen in the distributions P(q) of co-ordination 
numbers as depicted in Fig. 3. While P{q) is peaked at 
q = 6 for the Delaunay triangulations, it drops mono- 
tonically starting from q = 3 for the case of dynami- 
cal triangulations, and large co-ordination numbers are 
much more probable for the latter lattice type. It can be 
shown that for large co-ordination numbers the distribu- 
tion P{q) falls off as exp(— aging) with cr « 2 for Poisso- 
nian random lattices'*^, whereas for dynamical triangu- 
lations it declines much slower proportional to exp(— erg) 
with a = ln(4/3) « 0.3 (Ref. 50). 



III. THE RELEVANCE CRITERION 

Trying to decide for which kind of models disorder of 
the random-bond type constitutes a relevant perturba- 
tion, Harris* suggested the following line of reasoning. 
For a system of uncorrelated random bonds, the fluctu- 
ations (Jr{J) of the average coupling J in patches of the 
lattice decline with the linear patch size R according to 
the central limit theorem, that is, 



(5) 



where d denotes the spatial dimension of the system. 
Since disorder of the random-bond type couples to the 
local energy density, local fluctuations of the coupling 
induce such fluctuations of the effective local transition 
temperature, which decline identically with increasing 
patch size. Approaching the critical point t = {T — 
Tc)lTc = 0, the fluctuations in a correlation volume scale 
as 



a5(J)~e 



-d/2 ^ 



(6) 



where the power-law divergence ^ ~ t~'^ of the correla- 
tion length was used. For the critical behavior of the pure 
system to persist, these fluctuations should die out as the 
reduced temperature is linearly tuned to the critical point 
t = 0, i.e., one should have vd/2 > 1 or, with appeal to 
hyper-scaling, a < 0. In the converse case a > an al- 
tered universal behavior might me expected. Later on, 
Chayes et al}^'^'^ showed that under quite general condi- 
tions the specific-heat exponent of the disordered system, 
Q!dis, should be negative. This seemed plausible, since a 
was believed to coincide with the crossover or stability 
exponent of the respective renormalization-group (RG) 
fixed point^^'^^. Later on, however, it was claimed that 
in some cases the stability exponent might differ from a, 
such that even for positive a the regular critical behavior 
would prevail^°. Several examples of such behavior, and 
even the opposite case of disorder being a relevant per- 
turbation albeit a < 0, have been explicitly constructed 



for the case of hierarchical lattices^^"^'''. Thus, the gen- 
eral validity of the relevance criterion implied by Harris' 
argument has recently been the object of some debate, 
see, e.g., Refs. 66-68. 

As soon as the assumption of uncorrelated disorder de- 
grees of freedom is relaxed, Harris' reasoning is no longer 
applicable as it stands. For a random-bond model with 
long-range correlations of the bond variables, Weinrib 
and Halpcrin'^^ performed a rcnormalization group calcu- 
lation, which was later on refined^^. They find that, if the 
disorder degrees of freedom are correlated algebraically 
with a decay exponent a < d, disorder is irrelevant if 
a < 2 — Idja, whereas in the opposite case, the system 
flows towards a new, long-range correlated disorder fixed 
point. For more general disorder degrees of freedom and 
types of correlations including the ensembles of random 
graphs considered here, this argument can be adapted in 
the spirit of Luck's reasoning for the case of aperiodic 
structures'^ as follows. Consider a patch P of spherical 
shape with radius R and a volume of B{R) vertices*^ 
on a given realization of a triangulation. The average 
co-ordination number in P, 



J{R) 



B{R) 



(7) 



fluctuates around its expected value Jo = q [cf. Eq. (1)]. 

As the size of the patch is increased, R — > oo, these 
fluctuations decay as 



cjr{J) = {\J{R)-M)/Jo^{B{R)) 



-(l-<^) ^ ^-dh(l-w) 

defining the wandering exponent lo of the considered 
graph type. The Hausdorff dimension dh enters here to 
account for the cases of fractal graphs. In Eq. (8), the 
averages (•) are to be understood as the ensemble av- 
erages of the considered class of gra phs of a given total 
size. While for cj = 1/2 the usual l/ ^J{D{R)) behavior of 
uncorrelated random variables is recovered, for random 
lattices with long-range correlations of the co-ordination 
numbers one expects oj > 1/2, leading to a slowed-down 
decay of fluctuations. Near criticality, the fluctuation 
(T|(J) of the average co-ordination number in a correla- 
tion volume induces a local shift of the transition tem- 
perature proportional to \t\'^'^''^^~^^ ^IJ"^ . For the regular 
critical behavior to persist, these fluctuations should die 
out as the critical point t = is approached. This is the 
case when oj does not exceed the threshold value 



1 - 



a 



1 



dhV 2 — a' 



(9) 



provided that hyper-scaling is in effect. Conversely, for 
fluctuations satisfying uj > ujc a. new type of critical be- 
havior could occur. Re-writing Eq. (9) as 



1 - 2uj 



(10) 



it is obvious that for io = 1/2 the Harris criterion is 
recovered. 
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It is easily seen that the case of algebraicahy decaying 
correlations discussed in Ref. 31 is included in Eq. (10) 
as a special case. Specifically, for the case of connectivity 
disorder considered here, an isotropic power-law correla- 
tion would be given by 



Gqq{i,j) = {SqtSqj) - dist(i, j)" 



(11) 



where Sqi — — q denotes the co-ordination number 
defect at vertex i and the distance dist(i, j) is defined as 
the unique number of links in the shortest path of links 
connecting the two vertices i and j . Then, the fluctuation 
of the mean Squ over a spherical patch P of radius R is 
given by 



1 



(T\6q) 



const X R 



^ {5q,5qj) 
const X i?"" 



(12) 



as i? ^ oo. For short-range correlations a > dh the 
leading behavior is that of Eq. (5) , such that the Harris 
criterion stays in effect. For long-range correlations a < 
dh, on the other hand, the leading term is proportional 
to Comparing to Eq. (8), we arrive at the following 

expression. 



w = 1 — a/2dh, 



(13) 



such that from Eq. (10) we arrive at ac = 2 — 2dh/a, in 
agreement with the direct observation in Ref. 31. 

Since for systems with sufficiently long-range correla- 
tions of the disorder degrees of freedom oj > 1/2, such 
correlated disorder is more relevant than uncorrelated 
randomness in the sense that a change of universality 
class can already be expected for some range of negative 
values of a, cf. Eq. (10). If, on the other hand, corre- 
lations decay exponentially, the threshold ac = of the 
Harris criterion should stay in effect. 



IV. CORRELATORS AND WANDERING 
EXPONENTS 

In view of the discussion presented in Sec. HI, two 
complementary approaches towards a numerical determi- 
nation of the wandering exponent to present themselves, 
either a direct evaluation of the scaling of the fluctuations 
defined in Eq. (8) or an analysis of the correlation func- 
tion Gqq{i,j) of the disorder degrees of freedom defined 
in Eq. (11) to infer an estimate of u via Eq. (13). Both of 
these methods will be applied for the cases of Poissonian 
Voronoi'-Delaunay as well as dynamical triangulations. 
For both graph ensembles, a quenched average has to 
be taken over a number of graph realizations of the con- 
sidered ensemble. For this purpose, we generate lattices 
of spherical topology and consider the triangulations of 
varying co-ordination numbers as the basic objects and 
refer to their geometric duals as the "dual" lattices. 



Delaunay triangulations 
dynamical triangulations 
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FIG. 4: Comparison of the connected correlator Gqq{R) of 
the co-ordination number defects Sq = qi — q oi Poissonian 
Delaunay triangulations and dynamical triangulations of size 
= 125 000 triangles. The connecting lines are only drawn 
to guide the eye. The two short horizontal lines indicate the 
exact result /i2 = 10.5 (cf. Ref. 54) and the value Gqq{l) « 
— 1.2295 for the dynamical triangulations discussed in App. 
A, both valid in the limit A'^ — > oo. 



A. Analysis of Correlation Functions 



For the numerical determination of the correlator 
G'gq(i,j) a decomposition of the graphs into spherical 
shells around a given vertex is performed. This is done 
by first picking a vertex of the triangulation at ran- 
dom and a subsequent slicing of the graph into shells of 
equal geodesic link-distance around that vertex. For that 
given choice of initial vertex iq, two fundamental observ- 
ables can be measured, the volumes of the decomposition 
shells, 



(14) 



which, if properly averaged over, give the correlator of 
the unit operator, and. 



dist(2.'io) — -R 



(15) 



the average of which gives the correlator of Eq. (11). In 

view of Eq. (7) it is obvious that B{R) = Ef=o5ii('')- 
In the context of the dynamical triangulations approach 
to quantum gravity, there has been some debate on how 
to properly define (connected) correlators on an ensem- 
ble of random graphs, see, e.g., Ref. 69. The uncertainty 
concerns the order of taking the averages over a single 
graph and the graph ensemble, which has not been ex- 
plicitly specified in Eq. (11). The two possibilities are 
given by averaging the expressions of Eqs. (14) and (15) 
individually, i.e.. 



Gqq{R) 



(gy,(fi)) 
{9\\{R)y 



(16) 
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FIG. 5: Decay of the correlator —Gqq{R) of the co-ordination 
numbers of dynamical triangulations with TV — 500 000 trian- 
gles in a double logarithmic plot. The lines show fits of the 
functional form (18) to the data, where either the correction 
term was omitted {B — 0, thick dashed line) or included with 
B and b as free parameters (solid line). The ranges of the 
lines indicate the range of data points included in the fits. 

where the average (•) denotes a combined average over 
the starting vertices io and the graph realizations un- 
der consideration, which is obviously equivalent to the 
ensemble average of the graphs. The additional average 
over the starting vertices iq is merely included to improve 
the statistical accuracy in practical applications. On the 
other hand, one could also average on the level of the 
fraction, 

G„JR) = (^M^\. (17) 

It turns out that, at least for the quantum gravity graphs, 
these two ways of performing the average yield strikingly 
different results, even in the thermodynamic limit^^. The 
main motivation for using Eq. (16) is that an expression 
of the form {g^^^{R)), if integrated over all distances R, 
still yields some kind of susceptibility of the associated 
operator (f). Otherwise, however, correlators defined ac- 
cording to Eq. (16) behave rather pathologically. Defi- 
nition (17), on the other hand, corresponds to the nat- 
ural probabilistic definition of the average correlation of 
a given quantity at distance R, and is thus the unique 
"correct" definition in the given context and will be used 
throughout. Note, that in this case the average has to 
be performed carefully, since the maximum linear sepa- 
ration i?inax of two poiuts On the graph is not universal, 
but depends on the graph realization as well as on the 
chosen initial vertex. 

Figure 4 shows an overview of the short-distance be- 
havior of the correlator Gqq {R) defined according to Eq. 
(17) for Poissonian Delaunay as well as for dynamical tri- 
angulations as measured by averaging over 100 graph re- 
alizations sampled with one different starting vertex per 



FIG. 6: Finite-size scaling of the co-ordination number corre- 
lator ~Gqq[R*{N)] at distances R*{N) according to Eq. (21) 
for dynamical triangulations of sizes ranging from A'' — 1000 
to TV = 500 000 in a log-log plot. The solid line shows a fit of 
the functional form (20) to the data. 



1000 graph vertices each. The statistical errors here and 
in the following were determined via jackknifing (see, e.g., 
Refs. 70,71) over the 100 different graph realizations. Ob- 
viously, Ggg(O) — fj,2, corresponding to the values cited 
above in Sec. II. Note, that the negative correlation at 
distances R > reflects the fact, expressed in the Aboav- 
Weaire law (3), that a vertex with co-ordination num- 
ber q > 6 tends to have neighbors with g < 6 and vice 
versa. The correlation for nearest-neighbor points for the 
case of dynamical triangulations can be calculated in the 
thermodynamic limit by a series expansion of the results 
found in Ref. 54, yielding 0,5(1) « -1.2295, cf. App. 
A. This result is in reasonable agreement with the value 
found numerically for the example of TV = 125 000 trian- 
gles presented in Fig. 4, Gqq(l) « —1.34. The remaining 
difference gives a flrst indication of the presence of large 
finite-size corrections for the case of the highly fractal 
quantum gravity graphs. Going beyond i? = 1, a short 
glance at Fig. 4 reveals that the correlations are much 
more long-ranged for the case of the quantum gravity 
graphs than for the Voronoi-Delaunay random lattices. 
In fact, for the Delaunay triangulations we find that for 
the graph sizes up to TV = 500 000 triangles considered, 
co-ordination numbers of vertices at distances i? > 10 
apart are effectively uncorrelated up to the precision of 
our calculations. Due to this smallness of correlations it 
is not possible to determine their exact functional form 
from the accuracy of our measurements. However, a sim- 
ple exponential decay can be fitted reasonably well to the 
behavior found. Additional checks in support of this pic- 
ture will be presented below in this Section and in Sec. 
IV B. 

For the case of the quantum gravity graphs, we deter- 
mine the asymptotic behavior of the correlator (17) by 
fitting a suitably parametrized function to the numerical 
data. For the algebraic decay expected for this graph 
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TABLE I: Results of fits of the functional form (18) to the correlator Gqq{R) for dynamical triangulations with N = 500 000 
triangles and with various restrictions. 



Restriction 


A 


a 


B 


h 


-^min Rmax 


Q 


5 = 


-2.11(45) 


1.854(105) 







7 17 


0.981 


& = 3 


-2.37(29) 


1.938(92) 


5.6(19) X 10-=^ 


3 


4-37 


0.999 


6 = 4 


-2.04(21) 


1.841(75) 


1.27(40) X 10"® 


4 


4-37 


0.994 


none 


-2.30(24) 


1.919(76) 


3.1(67) X 10"^ 


3.15(57) 


4-37 


0.999 



type, we make the ansatz, 

Ggg{R) = AR-''{l + BR'>), (18) 

for i? > 1, taking into account an additional effective 
scaling correction with exponent b. From the discussion 
in App. B we expect 6 = = 4 to be a reasonable choice 
as long as no additional, more relevant, non-analytic cor- 
rections are present*^. Note that, since this "correction" 
term is more singular than the leading term 
this form is merely an effective description for a finite 
graph and distances small compared to the linear ex- 
tent of the graph. Figure 5 depicts the behavior of the 
correlator of co-ordination numbers for the case of quan- 
tum gravity graphs as well as fits of the form (18) to the 
data. We expect the range of applicability of the form 
(18) to be limited at both sides: for very small distances 
R, discretization effects are observed, whereas for large 
distances finite-size effects modify the expected behav- 
ior, i.e., higher-order terms neglected in Eq. (18). To 
account for these limitations, sampled points from both 
sides of the R range are successively dropped from the 
fit while monitoring the goodness-of-fit parameters 
resp. Q.^^ We find restricted fits with the correction ex- 
ponent b fixed at values 6 = 4 or, alternatively, 6 = 3, to 
match the data reasonably well. A full four-parameter 
fit with variable exponent b yields the intermediate value 
b = 3.15(57), in total indicating the presence of higher- 
order resp. non-analytic corrections. The fit results are 
compiled in Table I. Since we consider the correction 
term of Eq. (18) as an effective description, we take the 
fit with a variable value of the correction exponent b as 
the most reliable and quote as our best estimate of the 
decay exponent from this method 

a =1.919(76). (19) 

Due to the cross-correlations of Gqq{R) for different dis- 
tances R, the errors estimated by standard fit routines 
are biased; thus, instead, the errors were estimated by 
jackknifing over the whole fitting procedure. As a check 
of whether the used number of replicas in the disorder 
average is sufficient, we also performed the same analysis 
with subsets of the 100 realizations. Apart from the effect 
of rather sudden jumps of the fits into different, nearby 
minima, which always tend to occur with non-linear fit- 
ting procedures, we find all results completely consistent 



with each other within statistical errors. Thus, e.g., for 

the unrestricted fit of the form (18) we find (using the 
same ranges iimin--Rmax as before) a = 2.032(115) using 
only 50 graphs and a = 1.885(216) using only 10 graphs. 

A different way to determine the decay exponent a 
is based on a direct application of the finite-size scaling 
(FSS) behavior of the correlator. From FSS arguments 
it is known that the value of the correlator at a distance 
R* {N) scaled linearly with the size of the system behaves 
as 

Ggg[ii*(iV)] ~iV-"/'^\ (20) 

cf. App. B. The position of the reference points R* should 
be selected well within the scaling region of the consid- 
ered graph size. Here, we take the largest used graph size 
as a reference and choose 

where dh = 4. Obviously, the resulting distances R*{N) 
will in general be fractional numbers, for which no data 
are directly available. Since the accessible linear graph 
sizes N^/'^'^ are rather small, this discretization effect is 
quite pronounced. To circumvent this problem, we take 
a fit of the functional form (18) to the correlator with 
the correction term included as a suitable interpolation 
of the data used to estimate the value of Gqq[R* (N)]. 
As an estimate of the statistical error, we take the error 
computed for the closest integer distance. Figure 6 shows 
the resulting scaling plot and a fit of the functional form 
(20) to the data. We find the data not precise enough 
to resolve any corrections to the expected scaling form. 
The fit, including a series of graph sizes ranging from 
N = 1000 to A'' = 500 000, yields an estimate for the 
decay exponent of 

a = 2.09(13), (22) 

with quality Q = 0.53. For completeness, we have also 
performed the same analysis for the case of random Dc- 
launay triangulations. As expected, the resulting values 
of Gqq[R*{N)] do not show scaling according to Eq. (20), 
backing up our conjecture that the correlations are not 
algebraic in this case. 
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FIG. 7: Scaling collapse on the universal scaling function W 
defined in Eq. (B2) of the co-ordination number correlator 
Gqq (R) for dynamical triangulations of sizes A'^ = 20 000 to 
A'^ = 500 000 triangles. The vertical lines indicate the extent 
of the scaling window. 



Exploiting the FSS form (B2) of the correlator, it is 
also possible to re-scale the data such as to model the 
universal FSS function W introduced in Eq. (B2). Thus, 
plotting Gqq{R)N°-/'^'' as a function of the reduced dis- 
tance R/N^/'^'^ should, within the scaling region, produce 
data lying on a single master curve, irrespective of the 
lattice volume N under consideration. Such a scaling 
plot is shown in Fig. 7, where a nice scaling collapse of 
the data for an intermediate range of reduced distances 
R/N^/'^'^ is observed. For producing this scaling plot, 
we assumed a = 2, in agreement with the results found 
so far and further evidence to be presented below. The 
quality of the collapse does not change visibly on slightly 
changing the value of the exponent a to one of the values 
a — 1.919 or a = 2.09 found before. The extension of the 
scaling window in reduced distance R/N^^'^'^ can be very 
nicely read off from Fig. 7 to be 



0.1 < 



R 



^ 0.6, 



(23) 



yielding, e.g., 1 < i? < 7 for TV = 20 000 and 3 < i? < 16 
for N = 500 000, in perfect agreement with the observa- 
tions from the direct analysis of the correlator. 

Extending the analogy between the magnetic correla- 
tion function in a critical spin system and the correlator 
considered here, one notes that the decay exponent a 
defined above is related to the conventional critical ex- 
ponents hy a = 2 — dh — T]. A convenient and very precise 
method of determining ry is to consider the scaling of the 
integrated correlation function corresponding to the sus- 
ceptibility instead of analyzing the correlation function 
directly. The associated critical exponent j/i' then al- 
lows one to infer 77 via the scaling relation rj — 2 ~ "1 jv. 
Unfortunately, however, it follows from the relation (1) 
that the average co-ordination number on a closed trian- 
gulation is a fixed number, depending only on the number 




10 10 
1/<B(R)> 

FIG. 8: Decay of the averaged fluctuation (Jr{J) of co- 
ordination numbers against the inverse averaged integrated 
shell volume 1/{B{R)) according to Eqs. (7) and (8) for dy- 
namical triangulations (upper points) and Poissonian Delau- 
nay triangulations (lower points) . For both ensembles lattices 
with A'' = 500 000 triangles have been used. The solid lines 
show fits of the functional form (8) to the data. 



of triangles N. Thus the corresponding susceptibility al- 
ways vanishes identically, such that the method outlined 
above cannot be applied here. 



B. Averaged Fluctuations 

Instead of an analysis of the correlator of co-ordination 
numbers, the wandering exponent u> can be directly esti- 
mated by considering the scaling of the averaged fluctua- 
tions of co-ordination numbers and recurring to relation 
(8). This approach has the advantage of giving a nu- 
merical result also for the case of correlations decaying 
other than algebraically, as we saw for the Poissonian 
Voronoi'-Delaunay random lattices. We define the aver- 
age fluctuation J{R) as indicated in Eq. (7) with the 
average performed on the level of the fraction, in com- 
plete analogy with the definition (17) of the correlator. 
As before, the quenched average is performed over one 
starting vertex per 1000 vertices of the graph as well as 
over 100 different graph realizations for each graph size. 
Figure 8 shows the resulting decay of (J/j(J) against the 
inverse averaged integrated shell volume 1/{B{R)). The 
smaller number of data points for the case of dynamical 
triangulations results from the smaller effective linear ex- 
tent ^ N^/'^'^ due to their large fractal dimension. The 
relation of the two data sets nicely illustrates the much 
stronger correlations present in the dynamical triangu- 
lations model. From Eq. (8), we expect a linear decline 
of the curve in a logarithmic presentation, the slope be- 
ing given by 1 — w. On the other hand, corrections to 
the conjectured scaling behavior for very small distances 
R due to discretization effects as well as for large dis- 
tances comparable to the effective linear extent of the 
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FIG. 9: Finite-size scaling of the estimates of the wandering 
exponent u}{N) from fits of the functional form (8) to the 
data for dynamical triangulations of sizes A*' = 1000 up to 
N = 500 000 triangles. The line shows a fit of the form (25) 
to the data. 

lattices have to be taken into account. In Fig. 8 the scal- 
ing window of algebraic behavior according to Eq. (8) 
is nicely visible. The dramatic, exponential drop of the 
fluctuations as 1/{B{R)) approaches*^ l/{2 + N/2) is an 
effect of the topological constraint (1), as a consequence 
of which the fluctuation aii{J) vanishes identically at the 
maximum observed distance R = i?max- To obtain reli- 
able estimates for the wandering exponent co from a fit of 
the functional form (8) to the data, we again successively 
drop points from either side of the interval of distances R 
while monitoring the goodness-of-fit parameters resp. 
Q.^^ Note that, as before, due to the cross-correlations 
between the values of (Ju{J) for neighboring distances R, 
we have to resort to jackknifing over the whole fitting 
procedure to arrive at reliable error estimates. 

From the fits to the data for the maximum graph sizes 
N = 500 000 depicted in Fig. 8, we arrive at the following 
estimates for uj, 




0.50096(55), i? = 21, . . . , 41, Delaunay tr., 
0.72492(86), i? = 5, . . . , 14, dynamical tr.. 



(24) 

From the experience with the analysis of the correla- 
tor presented above in Sec. IV A, we expect rather pro- 
nounced finite-size effects still to be present at least for 
the case of the dynamical triangulations, such that the 
quoted statistical error certainly does not account for all 
of the deviation from the asymptotic result. Note that 
the result for Delaunay triangulations is perfectly consis- 
tent with a wandering exponent w = 1/2 resulting from 
either power-law correlations with an exponent a > 2, 
which, however, did not show up in the analysis of the 
correlator presented above, or an exponential decline of 
correlations in agreement with the observed non-scaling 
of the correlator. To account for the suspected addi- 
tional finite-size corrections present in the estimate (24) 
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TABLE II: Decay exponent a resp. wandering exponent u 
related by Eq. (13) for the ensemble of regular dynamical 
triangulations as estimated by various scaling methods. In 
each row, the directly measured value is printed in bold face. 



Method 


a 




Eq. (18) 


1.919(76) 


0.7601(95) 


Eq. (20) 


2.09(13) 


0.739(16) 


Eq. (25) 


2.021(78) 


0.7473(98) 


average 


1.987(76) 


0.7516(95) 



for dynamical triangulations, an additional FSS analysis 
is conducted by performing the analysis described above 
for the series of different graph sizes under consideration, 
ranging from N = 1000 to iV = 500 000. The result- 
ing FSS plot is depicted in Fig. 9. We expect a scahng 
approach of the following form, 

Lj{N) ^ ujoo + AN-'^ , (25) 

with an a priori unknown correction exponent 9. A fit 
of this form to the data yields: 

LOoo = 0.7473(98), 
A = -0.73(37), 

d = 0.264(70), ^ ' 

Q = 0.28, 

Note that the correction exponent is close to l/du = 1/4, 
which would correspond to an analytic scaling correc- 
tion. Fixing 9 ~ 1/4, we find the very similar result 
Woo — 0.7493(23). For the Poissonian Voronoi-Delaunay 
triangulations, on the other hand, we find only pretty 
small variations of the estimates U!{N) with the graph 
size, which are about of the same size as the statisti- 
cal fluctuations. Especially, for the largest two or three 
graph sizes, there is no visible drift between the esti- 
mates, such that we can safely take the result (24) for 
N = 500 000 triangles as our final estimate for the wan- 
dering exponent there. 

Finally, we also checked for the influence of the num- 
ber of disorder replicas on the results for the wander- 
ing exponent from the averaged fluctuations. For the 
Delaunay triangulations of size N = 500 000, using 50 
resp. 10 instead of 100 realizations yields estimates of 
uj = 0.50046(73) resp. uj = 0.49818(469), in very good 
agreement with the result for the full number of replicas. 
In the case of the dynamical triangulations graphs, reduc- 
ing the number of replicas to 50 resp. 10 and performing 
the fits (25) with the correction exponent 9 fixed at 1/4, 
we arrive at w = 0.7481(31) resp. uj = 0.7590(72). With 
an unconstrained correction exponent 0, we find the data 
with only 10 replica not precise enough to reliably apply 
the non-linear fit procedure. In Table II we collect the 
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final results for the decay exponent a resp. the wander- 
ing exponent lo for the case of dynamical triangulations 
from the various methods applied, using Eq. (13) with 
d;i = 4 to compute ui from the directly estimated decay 
exponent a or vice versa. 



V. CONCLUSIONS 

We have considered the applicability of a relevance cri- 
terion of the Harris-Luck type to the problem of coupling 
lattice models of statistical mechanics to the ensembles 
of Poissonian Voronoi-Delaunay random lattices and dy- 
namical triangulations resp. planar, "fat" 0^ Feynman 
diagrams. Following Luck's extension of Harris' origi- 
nal argument to the case of systems with aperiodicity, 
a relevance criterion is formulated for the case of ran- 
dom graphs with connectivity disorder, i.e., a random 
distribution of co-ordination numbers. Depending on a 
characteristic of the spatial correlations of the disorder 
degrees of freedom termed the wandering exponent uj, 
the threshold of relevance for a model with specific-heat 
exponent a predicted by the relevance criterion is shifted 
from Harris' value etc = to somewhere in the interval 
— cxD < ttc < 1, depending on the value of uj. 

To determine the values of the wandering exponent 
for the considered ensembles of random graphs, we have 
employed a detailed series of scaling studies. First, we di- 
rectly considered the behavior of the connected two-point 
correlation function of the co-ordination numbers of the 
vertices of the triangulations. For the Poissonian random 
lattices we find the correlations to decay very rapidly, and 
our analysis indicates that this decline is faster than any 
power of the distance. For dynamical triangulations, on 
the other hand, we find much stronger correlations, ex- 
hibiting a power-law decay within the scaling window. 
Due to the large fractal dimension of these lattices, how- 
ever, we find quite strong corrections to the leading scal- 
ing behavior to be present. Taking the correction terms 
carefully into account, we are able to determine the decay 
exponent a consistently with different methods, namely 
a direct analysis of the correlator, a finite-size scaling 
study of the correlation function at fixed distances and a 
scaling collapse of the data on a universal scaling func- 
tion. Via Eq. (13) this yields then an estimate for the 
wandering exponent lu. In addition to this analysis of 
the correlator, we investigate the scaling of the averaged 
fluctuations in spherical patches of increasing size, an 
expression which directly occurs in the definition of the 
wandering exponent. For the Voronoi'-Delaunay random 
lattices we find an estimate for the wandering exponent 
consistent with a value of w = 1/2 to high precision, 
which is in turn in agreement with the conjecture of cor- 
relations decaying faster than any power of the distance. 
Thus from Eq. (10), for Delaunay triangulations the pre- 
sented relevance criterion reduces to the Harris criterion 
ttc = 0, such that disorder of this type should be relevant 
for any model with positive specific- heat exponent. For 



dynamical triangulations, on the other hand, a combina- 
tion of the analyses presented allows us to conjecture^^ 
that the correlator decays with an exponent a = 2, lead- 
ing to a value of the wandering exponent of w = 3/4. In 
view of the relevance criterion (10), this leads to a rele- 
vance threshold as small as ac = — 2. As a consequence, 
the class of dynamical triangulations graphs should be 
a relevant perturbation for all known ordered models of 
statistical mechanics. 

There has been a number of numerical and analytical 
investigations of the effect of disorder from the classes 
of graphs considered here on the two-dimensional Potts 
model. For the case of a quenched ensemble of dynam- 
ical triangulations, a series of Monte Carlo simulations 
of the 5 = 2, 3, 4 Potts model with a = 0, 1/3, 2/3, 
respectively, showed a change of their critical exponents, 
indicating a shift to new universality classes'*'*'^^'^^. Ad- 
ditionally, it appears that the first-order transitions of 
the Potts models with g > 4 get softened to second-order 
ones as the model is coupled to a quenched ensemble of 
dynamical triangulations^^. For the case of percolation, 
corresponding to the — > 1 limit of the Potts model, for 
which due to its non-interacting character quenched and 
annealed averages coincide, an exact solution can be ob- 
tained by use of matrix model methods'"^. This solution 
as well as numerical simulations for this model''*, which 
has a — —2/3, also yield changed critical exponents, in 
agreement with the relevance threshold obtained here. 
For the case of Poissonian random lattices, simulations 
of the Ising or g = 2 Potts model, corresponding to the 
marginal case a = 0, yielded unchanged Onsager values 
for the critical exponents*^^"*^'^^; similar results where 
obtained for percolation^^, both results being in agree- 
ment with the relevance threshold ac = 0. On the other 
hand, for the g = 3, 4 Potts models with a > 0, Poisso- 
nian random lattices should be a relevant perturbation. 
However, an exploratory Monte Carlo study'''' as well as 
preliminary results from a high-precision series of Monte 
Carlo simulations of the authors''* for the three-states 
Potts model do not show any change of universal behav- 
ior, in contradiction with the relevance threshold pro- 
posed here. It remains an open question to be answered, 
e.g., by further high-precision analyses of the q = 3 case 
as well as the larger-a q = 4 Potts and Baxter-Wu''^ 
models, whether these findings are merely the effect of 
a crossover to new universal behavior occurring only for 
extremely large lattice sizes or whether there is possibly 
some physical reason for an argument of the Harris-Luck 
type not being applicable to the case of spin models cou- 
pled to Poissonian random lattices. For this purpose, a 
careful analysis of the counter-examples found in Refs. 
61-68 might be instructive. 
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APPENDIX A: NEAREST-NEIGHBOR 
CORRELATIONS FOR QUANTUM GRAVITY 
GRAPHS 

This Appendix is devoted to a short derivation of the 
result for the co-ordination number correlator at dis- 
tance one, Gqg(l), for the ensemble of quantum gravity 
graphs mentioned in Sec. IV A. Godreche et aL^^ con- 
sider topological correlations in the thermodynamic limit 
of the ensemble of (regular) dynamical triangulations via 
a generating-function technique. This allows them to 
compute the probability distribution Q/n of finding an 
edge connecting a vertex with co-ordination number I 
with a vertex with co-ordination number n by means of 
a series expansion. This quantity is related to the prob- 
ability of finding an Z-vertex in the neighborhood of an 
n-vertex, which is of interest here, as follows: 



where 



p, -l^n 

Z In 



(Al) 



(A2) 



Then, the distance- 1 correlator can be expressed as 
G,,(l) = (6-0(6 -n) = ^{1-12S + Z), (A3) 



the symmetry property Qin = Qnl) yields Z w 0.96697. 
Hence, we arrive at 



Ggg{l) = 36 



Z-1 



-1.2295. 



(A6) 



APPENDIX B: FINITE-SIZE SCALING OF THE 
CORRELATOR 

In this Appendix we give a short justification of the 
finite-size scaling (FSS) method used in Sec. IV A to de- 
termine the scaling exponent of the two-point correlator 
of co-ordination numbers and the correction terms taken 
into account. We consider making a real-space renor- 
malization transformation with a re-scaling factor b (see, 
e.g., Refs. 1,80,81). After n iterations of the re-scaling, 
the two-point correlation function of some operator un- 
der consideration can be written as 



G^^(i?;t,^) = 6-^-*"G^ 



6" N 



(Bl) 



where denotes the scaling dimension of the operator 

(j) and yi the temperature-related renormalization-group 
eigenvalue. Stopping the scale transformation at an it- 
eration such that N-'^b'"^'^ = K, i.e., 6" = {KNy/^'^ = 
{N/Noy/'^'', we arrive at 



N 



W 



R 



iN/No)^_ 



(B2) 



introducing a universal scaling function W. Therefore, if 
R is scaled linearly with N^/'^h^ e.g., R* = N^/^^^ /2, the 
correlation function scales as 



(B3) 



where S = J2in Qin/l- Integrating the generating func- 
tion of Qln, 



Q{x,y) = ^QinX^y", 



one arrives at: 



JO X 7 „ ^ 



X 

Z = 36 (\y^M_m 

Jo y 



36 E 



Qln 

In 



(A4) 



(A5) 



The first integral can be performed exactly, yielding S = 
P(l) = 1/6. The double integral in the second line of Eq. 
(A5) could not be evaluated in closed form. Instead, the 
series expansion of Q(a;, y) used in Ref. 54, performed up 
to orders I < 50 and n < 100 (and additionally exploiting 



Now, instead of doing FSS, consider the finite-size 
"field" as a scaling correction and stop the re-scaling at 
i?/6" = K, 



G^^{R;t, jf) — 



R 

i?0 



F 



R_ 

Rq 



l_ fR_\ ' 
'N\R^) 



where F is another scaling function. Hence, at criticality 
one has 



G^^{R; 0, 



R_ 
Ro 



N \R^) 



(B5) 

which should be a reasonable approximation as long as 
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All distances on the graphs considered in this paper are to 

be understood as the unique number of links in the shortest 

path of links connecting two vertices. 

An implicit assumption in applying such scaling arguments 
is that the ensemble of considered random graphs behaves 
as a system of statistical mechanics at a critical point. From 
the dynamical triangulations approach to quantum gravity 
it is known that this is indeed the case, see, e.g., Ref. 38 
Note that due to the correlations between values of Gqq [R) 
for different distances R, the absolute values of resp. Q 
are not immediately meaningful; relative changes, however, 
are. 

*® The number of vertices of a spherical triangulation consist- 
ing of N triangles is iVo = 2 + N/2. 

*® Note that all critical exponents of the dynamical triangu- 
lations model, which could be calculated exactly, where 
found to be rational numbers, see, e.g., Ref. 38. 



